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Abstract. We investigate the dust composition of detached shells around carbon stars, with a focus to understand the origin of 
the cool magnesium-sulfide (MgS) material around warm carbon stars, which has been detected around several of these objects 
iHonv et al.l2002li . We build a radiative transfer model of a carbon star surrounded by an expanding detached shell of dust. The 
shell contains amorphous carbon grains and MgS grains. We find that a small fraction of MgS dust (2% of the dust mass) can 
give a significant contribution to the IRAS 25 /im flux. However, the presence of MgS in the detached shell cannot be inferred 
from the IRAS broadband photometry alone but requires infrared spectroscopy. 

We apply the model to the detached- shell sources R Scl and U Cam, both exhibiting a cool MgS feature in their ISO/SWS 
spectr a. We use the shell p arameters derived for the molecular shell, using the CO submillimetre maps iLindavist et alll999l: 
ISchoier & Olofssorl200ll) . The models, with MgS grains located in the detached shell, explain the MgS grain temperature, as 
derived from their ISO spectra, very well. This demonstrates that the MgS grains are located at the distance of the detached 
shell, which is a direct indication that these shells originate from a time when the stellar photosphere was already carbon-rich. 
In the case of R Scl, the IRAS photometry is simultaneously explained by the single shell model. In the case of U Cam, the 
IRAS photometry is under predicted, pointing to a contribution from cooler dust located even further away from the star than 
the molecular shell. 

We present a simple diagnostic to constrain the distance of the shell using the profile of the MgS emission feature. The emission 
feature shifts to longer wavelength with decreasing grain temperature. One can therefore infer a temperature and a correspond- 
ing distance to the star from the observed profile. Such a diagnostic might prove useful for future studies of such systems with 
SIRTF or SOFIA. 

Key words. Stars: AGB and post-AGB - Stars: carbon - Stars: mass-loss - Stars: individual: R Scl & U Cam- Circumstellar 
matter - Infrared: stars 



5^ ! 1. Introduction ments jLattimer et al.lll978l ILodders & Feglevlll999l) . but are 

chem ically unstable in oxidising, i.e. oxygen-rich, surround- 
The so-called "30" ^m emission feature is commonly de- ings J Num et alJllQsi . The featur e profile varies con sider- 
tected in the infrared (IR) spectra of carbon-rich (C-rich) ablyjrom source t o source (e.g . iGoebel & Moselevl f l985 } 
evolved stars, from low to intermediate mass-loss asymptotic | Waters et a1 J iHrivnak et all l200(t . iHonv e~tm& 



giant branch stars (AGB) through infrared carbon stars and haye presented a detailed spectral comparison between the 
post-AGB stars to planetary nebulae (e.g. |Forrest et al. [l98ll emission expected from MgS grains and the profiles of the 



Omontl ll993l lYamamura et alJ Il998l Hony et alJ l2002h . The 



"30" fim emission feature found in the IR spectra of C-rich 

"30" urn feature, which exhibits itself as a prominent broad evolved stars obtained wi th the Short Wavelength Spectrograph 

emission excess in the IR spectrum, extending from -24 to (sws ^ | de Graauw et all fl99(3) on b oard the Infrared Space 

-45 //m, is commonly attributed to magnesium-sulfide (MgS) observatory (ISO. lKessler et alll99fil) . Although there are also 

dust grains (e.g. | Goebel & Moseley l [1985J | Begemann et alJ some profile var i ations tha t can be attributed to grain-shape 

H994J). MgS grains are expected to form in C-rich environ- variationSi it was demonstrated that the main variations in the 

appearance of the "30" fim feature between distinct sources can 

Send offprint requests to: S. Hony, e-mail: shony@rssd.esa.int be understood as a result of differences in the (average) MgS 

* based on observations obtained with ISO, an ESA project with in- grain temperature. Conversely, one can derive the average grain 

struments funded by ESA member states (especially the PI countries: tem pe rat ure from the observed profile. One of the interesting 

France, Germany, the Netherlands and the United Kingdom) with the r. £ lU . , j t u j- c i «i_ ■ 

... „i' . _ ,» T .„. findings of that study was the discovery of several hot , i.e. 
participation of ISAS and NASA & 3 J 
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Fig. 1. The ISO/SWS spectrum of R Scl. An example of an op- 
tically bright carbon star that exhibits the "30" //m feature at- 
tributed to cool MgS grains. The "30" /mi feature is shaded in 
grey. At shorter wavelengths there are several prominent pho- 
tospheric, molecular absorption bands. The emission band at 
1 1.3 /im is attributed to silicon-carbide grains. 



optically bright, carbon stars that exhibit a cool MgS profile. 
In Fig. [l] we show the spectrum of R Scl. R Scl is an example 
of such an optically bright carbon star with a very prominent 
"30" fim feature. One of the most likely explanations for this 
phenomenon is the presence of MgS grains far from the star, 
whereas there is a lack of MgS grains closer by. 

Some optically bright carbon stars are known to possess a 
detached shell of gas and dust a round them fe.g. lWaters et alJ 
Il994l l Schoier & Olofsson 2001). In this paper, we explore the 
effect that MgS grains located in such a detached shell will have 
on the IR spectra and photometry of these stars. We build a 
grid of radiative transfer models of a carbon star surrounded by 
a detached dust shell which contains amorphous carbon (a-C) 
and MgS grains. We compare the synthetic spectra from these 
models with the IRAS photometry of known carbon stars. We 
further apply the model to two specific carbon stars: R Scl and 
U Cam. The detached shells aroun d these sources are clearly 
resolved in molecular line emission ( Schoier & Olofsson 2001 ) 
and good quality ISO/SWS spectra are available to compare 
with the synthetic spectra. This allows us to address the ques- 
tion on whether the MgS grains are located at the same distance 
as these molecular shells. 

1.1. The origin of the dust far from the star 

Carbon stars are known to loose large amounts of envelope 
material through a dust driven wind. The observational indi- 
cations that this mass loss may occur in a non-steady fash- 
ion are important for the present study. Perhaps the most 
complete picture arises from the IRAS colour-colour diagram 




[12]-[25] 



Fig. 2. IRAS colour-colour diagram of galactic carbon stars (di- 
amonds). Stars that exhibit a single temperature SED are shown 
in black. Stars with a far-IR excess are grey. We also show the 
tracks that the models follow in the diagram. The tracks are 
for a mass-loss burst lasting for 200 years. The dust mass-loss 
rate increases from the bottom to the top curves (2 10~ 8 -2 10~ 7 
Mo/yr). The time elapsed since the burst is indicated in the fig- 
ure. The difference between the dashed and the full lines is the 
relative amount of MgS. At the onset of the burst, when the 
dust is warmest, the contribution of MgS makes the tracks dif- 
fer substantially. 

(Ivan der Veen & Habingll988l) of C-rich evolved stars (Fig. [2). 
This diagram shows the 60 yum over 25 /im flux-ratio (expressed 
as a magnitudes difference) versus the 25 /mi over 12 /mi flux- 
ratio based on the IRAS photometry. A star's position in this 
diagram is thus a measure of the shape of its IR spectral energy 
distribution (SED). The diagram shows that most of the carbon- 
star SEDs are dominated by a single temperature. However, 
about one third of them show evidence of a second, cooler 
dust component. This cooler dust component manifests itself 
through a far-IR excess. 

The common interpretation is that this excess is due 
to a detached dust shel l located away fr om the star (e.g. 
IWillems & de Jon""lll988t Izuc kermanl 1 1 99 3h . This shell orig- 
inates from a previous epoch of mass loss, when the mass- 
loss rate exceeded the present day mass-loss rate by several 
factors. Although there has been some discussion in the lit- 
erature on whet her this far-IR excess is actually relat ed to 
the carbon stars dlvezic & ElitzuJl 9*951 lEgan et all 19 96). sev- 
eral of these detached shells have now been direc tly detected. 
These shells have been re s olved either in the IR faoungetalj 
[ 19931: IWaters et alJ Il994l llzumiura et all Il995i 1 1 99(1 1 19971: 
llzumiura & Hashim oto 2000) o r in molecular lines in the 
submillimetre domain lYamamu ra et alj"l993l lOlofsson et alJ 
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119961: iLindavist et aljfl99^ ISchoier & Olofssonl Eool . In a 
few cases, the optical light scattered by the material in 
these shells has been directly detected Jlzumiura et alJ feoOO: 
IGonzalez Dekado et all200ll) . 

The origin of these shells, i.e. the cause of the change in 
mass-loss rate, is matter for debate. The earliest suggestion was 
that occurrence of these shell was closely related to the photo- 
spheric composition of the AGB star. Specifically, it was sug- 
gested that these shells arise as a result of the thermal pulse, that 
increased the number ratio of C-atoms to O-atoms in the pho- 
tosphere of the AGB star from less than unity (O-rich) to la rger 
than u nity (C-rich) JWillems & de JongHl988t IChan & Kwokl 
Il988l) . Later studies showed that these shells may arise due to 
thermal pulses in general, as long as the star is close to the 
tip of the AGB, i.e. with a luminosity close to the critical lumi - 
nosity (e.g. lSchroder et al.ll998l:ISteffen & Schonberneil200Cl) . 
Recent numerical simulations of the momentum transfer in the 
dust driven winds of AGB stars show that modulations in the 
frictional coupling between the gas and the dust can also caus e 
non-steady mass loss and shell formation ( Si mis etal1l200ll) . 
These results have been obtained for a C-rich gas and dust com- 
position and it is at present not clear whether the same mecha- 
nism applies to O-rich environments. 

1.2. The composition of the material in the shell 

Closely related to the question of the origin of these detached 
shells is the question of their composition. The molecular and 
dust composition in the shells is a tracer of the chemical and 
physical conditions, in particular the C/O number ratio, in 
the atmosphere of the star at the time the shell was formed. 
Ther e are various indirect eviden ces that these shells are C-rich. 
Buiarrab al & Cernicharollfl9 94) find that the molecular line ra- 
tios for the various molecules are indicative of a C-rich chem- 
istry. From modelling of the SED, one also tends to find that the 
opacities of C-rich dust agree better with the obser ved far-IR 
excesses than O-rich dust (e.g. lBagnulo et al1ll997l) . However 
these models are not conclusive and models with O-rich dust 
can also be fitted to the SEDs. The possibility of detecting MgS 
grains in these shells is a tantalising one, since this gives us a 
direct handle on the prevailing C/O ratio at the time the matter 
was ejected. This is due to the fact that the "30" /mi feature, 
attributed to MgS grains, is only found in C-rich environments 
and has, until now, not been detected in O-rich environments. 

The paper is organised as follows. In Sect.[5]we present the 
description of the model. The results of the model are presented 
and compared to IRAS observations in Sect. [3] We apply the 
model to R Sculptoris and U Camelopardalis in Sect. 0] Both 
these sources have a "30" /im feature in their ISO/SWS spectra 
and the detached gas shells have been spatially resolved with 
molecular line observations. This allows us to test whether the 
molecules and dust are co-spatial. Finally, in Sect. [5] we sum- 
marise our findings and discuss future prospects. 

2. Model description 

Our model consist of a single spherical dust shell surrounding 
a carbon star with a low current day mass-loss rate. The radi- 



ation of the star is partially absorbed by the dust grains in the 
shell. The absorbed light heats the dust grains, which will sub- 
sequently emit light at IR wavelengths. 

The optically visible C-star is represented by a Planck func- 
tion of 2500 K. This temperature is representative for the en- 
semble of optically bright c arbon stars with co ld MgS grains 
in the ISO/SWS sample from lHonv et al.N2002l). where w e use 
the effective temperatures as given by Berg eat et alJfcOOlh The 
luminosity of the star is set to 10 000 L , which corresponds to 
a stellar radius of 535 R Q . Of course, the whole system can 
simply be scaled to a different luminosity for the central star by 
scaling all dimensions with the appropriate value. 

2.1. Shell parameters 

We assume a constant expansion velocity, with a typical value 
of 15 km/s. The expansion velocity is constant. During the burst 
that produced the shell, the mass-loss rate was constant. Before 
and after the burst the mass loss was negligible. Due to the 
expansion, the dimensions of the shell are related to the time 
since termination of the burst as: 

R in = /?* + v exp * t and R out = R in + v exp * At, (1) 

where is the radius of the star, /?;„ and R out are the inner and 
outer edges of the dust shell after time t . At is the duration of the 
high mass loss. We calculate synthetic spectra for times after 
ejection (f) from 300 to 10 000 year. For t < 300 year the opti- 
cal depth through the dust shell is very large and these systems 
will look essentially like high mass-loss carbon stars, instead 
of detached-shell sources. For t > 10000 year the MgS grains 
are too cold to have any detectable signature in the resulting 
SED. We vary Af between 200 to 1500 year. The corresponding 
values for the shell thickness (1 10 16 - 7 10 16 cm) bracket the 
thickn ess of the detached shells derived in Sch oier & Olofssonl 
(12001 1) 

The gas density in the shell follows from Pf, ds {R) = 
M/(47r/? 2 v eX p), where p gas is the gas density as a function of 
radius (R) and M is the mass-loss rate. The dust density (pdust) 
is important for the radiative transfer calculations. This value 
is the ratio of the gas density to the gas-to-dust mass ratio Q¥), 
which implies that the models are determined by the dust mass- 
loss rate (Mdust = M/W). We calculate model spectra for the 
dust mass-loss rate in the range 2 10~ 9 - 2 10~ 7 M Q /yr. This 
translates to a total mass loss-rate in the range 5 10~ 7 - 10~ 4 
Mo/yr, using a typical value of V F=250. However, it is not clear, 
whether such a standard value is always applicable (see also 
Sect.g). 

2.2. Dust parameters 

We use a-C and MgS as the compositi on of the dust. The op - 
tical properties of a-C are taken from Preibisc h et alJ (11993). 
We use an a-C grain-size distribution according to: n(a) oc 
a -3 - 5 with 0.01 pmn < a < 1.0 fan, where n{a) is the number 
density of grains with radius a. For the MgS grains we use the 
same optical properties used to fit the observati ons of the C-rich 
post-AGB star HD 56126 faonv etalJl2003h . These proper- 
ties are based on the calculated absorption cross-sections in the 
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IR, u sing the optical constants as published by Bege maiin et alJ 
(1994) We calculate the IR absorption cross-sections for a con- 
tinuous distribution of ellipsoids shape-distribution. This dis- 
tribution results in a resonance which agree s best with the 
observed feature (e.g. iBegemann et aD Il994l ISzczerba et alJ 
19991 iHonv et alJl2002h 7 Unfortunately, the optical properties 
of MgS have not been measured below 10 fim. We assume a 
constant absorption cross-section below 1 /jm of 7rxl0~ 4 fim 2 
for a grain radius of 0.01 fim. We further assume a linear de- 
crease in the absorption cross-section to from 1 to 2 fim and 
between 2 and 10 fim. The cross-section assumed here is identi- 
cal to the one used to model HD 56126. The location of the dust 
around that sta r is well determin ed from extensive IR imag- 
ing st udies (e.g. lDaval et alll998lljura et alfeOOollKwok et alJ 
2002) as the dust shell is clearly resolved. The assumed absorp- 
tion cross-section below 1 fim yields a MgS grain temperature, 
at the location of the dust shell of HD 56126, which agrees well 
with its SWS spectrum faonv et all2003l) . 

We vary the relative proportion of MgS to a-C from 2 to 
10 percent by mass. A smaller fraction of MgS will not pro- 
duce a significant "30" /urn feature. A larger fraction of MgS 
is unlikely to be present because the amount of MgS grains 
is limited by the abundance of Mg and S atoms. For a solar 
composition gas, the maximum mass contained in Mg+S pairs 
relative to the mass of the C-atoms is 20 percent. However, not 
all available Mg, S and C atoms are necessarily be condensed 
into dust grains. 



2.3. Dust radiative transfer code 



We use the proprietary dust radiative transfer code MODUST. 
This code solves the monochromatic radiative transfer equa- 
tion, from UV/optical to millimetre wavelength, in spherical 
geometry subject to the constraint of radiative equilib rium, us- 
ing a Feautrier type solution method (Feautrier 1964; Mihalas 
1978). This yields the temperature of the dust grains. The code 
allows to have several different dust compon ents of various 
grai n sizes and shape s. We refer to iBouwman et alJ (|2000) 
and |Bouwman| ( 1200 II) for a description of techniques used 
in MODUST. MODUST yields both the model SED and in- 
tensity maps. From the obtained model spectra and intensity 
maps, it is straightforward to simulate observations. For ex- 
ample, we simulate the synthetic SWS spectra by overlay- 
ing the apertures as a f unction of SWS subband, as listed in 
Ide Graauw et alJ dl996l) . on the intensity maps. IRAS broad- 
band fluxes are simulat ed by convolvin g the spectra with the 
IRAS bandpasses (Beichman et alJl988l) . We solve the transfer 
equations at 36 wavelengths, on a logarithmically spaced wave- 
length grid, from 0.2 to 1000 fim. This is sufficient for correctly 
integrating over the wavelength, when performing the radiative 
equilibrium calculation. These wavelengths also cover the en- 
tire SED. In addition, we calculate the synthetic mid-IR spec- 
trum from 2 to 45 fim, with linear wavelength steps of 0.5 /mi, 
to compare with the observed SWS spectra. 




1 10 100 1000 

f 12 [Jy] 

Fig. 4. IRAS colour-magnitude diagram of galactic carbon 
stars. The thick black line shows the mean of the [25] -[60] 
colour as a function of 12 fim flux for the carbon stars (shown 
in grey). The sources with the largest 60 fim excess are all very 
faint at 12 fim (to the left of the dashed line). The photometry 
of these sources is probably affected by cirrus contamination or 
confusion due to the poor angular resolution of IRAS. 

3. Results 

In Fig. we show a subset of the SEDs that result from our 
model. We show the main effects of varying the mass-loss rate, 
the fraction of MgS grains and the duration of the burst. We 
also indicate which part of the SED contributes to the various 
IRAS photometric filters. There are a few general points that 
emerge from the analysis. 

1. The "30" fim feature is detectable even when present at 
a concentration of only 2 percent. This is especially true 
shortly after the ejection, when the MgS grains are rela- 
tively warm and the spectra clearly exhibit the characteris- 
tic onset at ~26 fim. 

2. MgS contributes considerably to the IRAS 25 fim flux dur- 
ing the first 2000 year when it is warm enough to emit 
within the IRAS 25 fim bandpass. 

3. The shell becomes quickly optically thin for the radiation 
of the central star and therefore the emission of MgS scales 
with the amount of MgS in the shell. 

4. The 12 fim excess disappears within ~750 year. 

5. As expected, the IR excess increases with larger shell mass 
and is more pronounced for a short fierce burst than for a 
prolonged ejection at a lower mass-loss rate, even for the 
same amount of total material shed. 

3.1. The IRAS colour-colour diagram 

In Fig. |2 we compare the tracks as predicted from our 
model with the IRAS colours of known galactic carbon 
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dust mass-loss rate=2 10 
MgS=2% 
At burst =200 




X [jim] 



Fig. 3. Synthetic SEDs from the expanding dust shell containing amorphous carbon and MgS grains. In each panel, the curves 
from top to bottom represent the system 300, 600, 900, 1200, 1500, 2000, 3000, 5000 and 10 000 year after ejection of the shell. 
The top panels show the effect of varying the percentage of MgS grains. In the bottom six panels, the mass-loss rate increases 
from top to bottom and the duration of the burst increases from left to right. Also indicated are the wavelength regions that 
contribute to the different IRAS broadband filters. The MgS grains are the cause of the "30" /J.m feature, the broad extra emission 
band between ~23 and ~45 //m. 
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Fig. 5. Illustration of the definition of the centroid position of 
the "30" fim feature. We shows two model spectra of a detached 
shell containing and 10 percent MgS. The difference is the 
"30" fim feature due to the presence of the MgS grains. On 
this scale equal surface corresponds to equal energy. The cen- 
troid position is indicated by the dashed line, which divides the 
"30" fim feature into equal surface. 



stars. The diagram is constructed in the following way. We 
cross-correlate the positions of th e sources in the G eneral 
Catalog of Galactic Carbon Stars JAlksnis et al.ll200ll) with 
the sources within 30" in the IRAS poin t source catalogue 
Jjoint IRAS Science Working Groudll988l) . This yields 3417 
associations. From these sources we only use those with a rela- 
tive flux uncertainty less than 50 percent in the IRAS 12, 25 and 
60 fim filter. In Fig.|2] the maximum uncertainties are 0.55 and 
0.6 and the average uncertainties are 0. 1 1 and 0. 15 in [12]-[25] 
and [25] -[60], respectively. 

The tracks explain the general features of the observed 
colour-colour diagram well. The colours of the systems imme- 
diately after ejection of the shell are similar to the colours of 
the infrared carbon stars, because of the large optical depth in 
the shell. The expansion of the shell causes the optical depth to 
rapidly decrease and the 12 fim excess to diminish as the dust 
becomes cooler. Therefore, the source moves rapidly to the left 
in the diagram. The stars spend by far the longest time on the 
vertical track, at [12]-[25] « -1.2, while the 60 fim excess due 
to the cool dust slowly diminishes. This explains the cluster- 
ing of the stars in this region of the diagram. The effect of the 
"30" fim feature due to MgS grains can only be noticed within 
2000 years after ejection. However, the influence of MgS on the 
colours is not unique, because the effect of the MgS grains in 
the colour-colour diagram can be mimicked by a higher mass- 
loss rate or a longer burst duration. This implies that, due to 
the presence of MgS grains in these shells, and thus of the "30" 
fim feature, the mass in the shell may be overestimated, when 
determined from the IRAS photometry alone. 

We note that there are several carbon stars in the IRAS 
database with 60 fim excesses larger than the ones shown in 
Fig. 13 We have not attempted to reproduce these IRAS colours 
with our model. It is likely that there is a sizeable contamina- 
tion of the far-IR fluxes among t he carbon stars, as has been 
argued by llvezic & ElitzuJ Jl995t) . Specifically, we point out 
that the largest IRAS excesses are found for the sources with 
low 12 fim fluxes, i.e. sources far away with a much larger 




1000 10000 
shell distance [R.] 

Fig. 6. MgS temperature and the centroid position of the 
"30" fim feature as a function of distance to the star. The di- 
amonds show the MgS temperature at the inner edges of the 
detached dust shell. The squares represent the centroid position 
of the "30" fim emission feature versus the typical distance of 
the shell, which is defined as the inner radius plus one quarter 
of the shell thickness. The grey dashed line shows a simple an- 
alytic representation of the relation between the distance and 
the centroid position (see Eq. (0). 



chance for confusion or cirrus contamination. The carbon stars 
with IRAS colour [25] -[60] > are almost exclusively found 
among sources with F12 < 5 Jy (see Fig. |4). In this paper, we 
are only concerned with the sources that have [25]-[60] < 0.5 
and clearly have a detached shell associated with them. In or- 
der to obtain 60 fim excesses in the range [25]-[60] - 1 or even 
higher (see Fig.@J, we have to increase the total dust mass to 
1O~ 4 M , which in turn yields an estimated mass-loss rate an 
order of magnitude higher than those derived for the best stud- 
ied cases. It is not clear whether the same detached-shell model 
can be invoked to explain such extreme excesses. 

3.2. MgS as a temperature/distance diagnostic 

Because the shell becomes rapidly optically thin, we can give 
an approximate temperature of the MgS grains as a function of 
shell radius and effective temperature of the star. The effects 
of the expansion of the shell are shown in Fig. |6] We show 
two curves. The diamonds show the dependence of temperature 
on distance. The squares show how this temperature translates 
into a shifting of the "30" fim feature (see also Fig. The 
shift is expressed in the centroid position of the MgS feature, 
i.e. the wavelength which divides the feature in equal energy 
halves, see also Fig.|5] This quantity is in general more easily 
derived from the spectra and is also a better diagnostic of the 
temperature than the peak position. For the distance of the shell 
we take a "typical" distance, defined as the inner radius plus 
one quarter of the shell thickness. This is to account for the 
fact that the dust on the inside of the shell is hotter and thus 
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Table L Shel l parameters of R Scl and U Cam. "Taken from 
ISchoier & Olofssoi] J200lh . *The dust parameters for U Cam 
only concern the dust located inside the SWS aperture (see text 
for details). 



.''IRAS 60 



IRAS 12,25 




brighter than on the outside. The points shown in Fig.[6]are for 
shell thickness between 250 and 2000 which corresponds 
to a burst duration from 200 to 1500 years. An approximate 
power-law expression (valid for 30 < Aq 30 < 36 yum) is given 
by: 



3D 



R, 



y2l jjon 



\25WK) 



(2) 



where /? s h e ii is the typical distance of the shell. Note that, in 
physical terms, the shape of the emission feature is a weighted 
sum of Planck functions, depending on the MgS grain temper- 
atures, convolved with the wavelength dependent MgS absorp- 
tion cross-section. Therefore, the simple power-law expression 
and the constants in the first term of Eq. (0 lack a direct phys- 
ical interpretation. With optical photometry of the central star, 
both T+ and /? s h e ii (in arcsec) can be determined independent 
of distance, as the angular size of the star (/?* in arcsec) is set 
by and the photometry alone. How practical Eq. will 
be in reality remains to be demonstrated, because it is known 
that there are also some differences in the shape of "30" fim 
feature between sources that are not caused by temp erature but 
are m ore likely due to the grain-shape distribution ( Honv et all 



2002). Such a difference can shift /lc,3o by up to 2.5 fim: a fac- 
tor 4 in shell radius! Because the main deviation in the profile 
is a 26 fim excess, the tendency will be to find a lower value of 
Ac,30 and therefore derive a too small shell radius. 

4. Application to sources with ISO/SWS spectra 

In the previous section we have shown that the presence of 
MgS in the detached shells can have a significant effect on the 
IRAS 25 /mi photometry. However, the presence of MgS grains 
in the shell cannot be inferred from the IRAS photometry alone. 
In this section we focus on those warm carbon stars that show 
a "30" fim feature in their SWS spectrum. 



4.1. Application to Ft Sculptoris 

We first apply the model to R Scl. This star is among the best 
studied optically bright carbon stars with a far-IR excess. Its 
detached shell has be en resolved in molecular lines emission 
(Olofsso n et al.lll996l) and scattered light from the detached 
shells is detected dGonzalez Delgado et al .2001) . Furthermore, 
R Scl has the strongest "30" fim feature of the optically bright 
carbon stars observed with SWS. 



Fig. 8. The angular dimensions of the circumstellar shell of 
R Scl compared with the regions probed by ISO/SWS and 
IRAS. The black rectangles indicate the dimensions of the two 
most relevant SWS apertures. The dashed circles give an indi- 
cation of the regions corresponding to IRAS 12, 25 and 60 fim 
point source measurements. The circle for the IRAS 100 fim 
filter lies outside the figure. The shell extends from 13" to 
26", which corresponds to 7 10 16 - 1410 16 cm at a distance 
of 360 pc. The complete shell is probed by the IRAS measure- 
ments. The SWS instrument detects only a part of the shell ra- 
diation. Of course, the back and front side of the spherical shell 
do contribute to the SWS flux. 



For the parameters of the dust shell (see Tabled, we take 
the inner and outer edge of the shell o f the detached gas shell as 
listed by Sch oier & Olofssonl d200 lL their Table 6). These pa- 
rameters are derived from mapping observations of the molecu- 
lar CO emission, in which the shell is clearly resolved. We use 
the same distance to the star as these authors do, to facilitate 
a direct comparison of the parameters, which they derive and 
ours. With the inner and outer edge of the shell fixed, we have 
only two free parameters: the dust mass-loss rate and the mass 
fraction of MgS. We vary these parameters and compare the re- 
sults with the IRAS and ISO observations. The strength of the 
IR continuum excess is set by the amount of amorphous car- 
bon in the shell, which scales linearly with the dust mass-loss 
rate. The fraction of MgS influences the strength of the resul- 
tant "30"yimi feature. The best agreement with the observations 
is found with 4 percent of MgS and a dust mass-loss rate of 
only 8 lO~ 9 M /year. This best fit is shown in Fig.0 We de- 
rive a total dust mass in the s hell, of 1 .2 10~ 5 M B . The gas mass 
as determined by ISchoier & Olofssonl fcOOlh is 5.4 10"" 3 M Q , 
which yields for the gas-to-dust mass ratio: *F=450. The model 
parameters are given in Tabled 

We show a sketch of the model, with the dimensions of the 
relevant SWS apertures and the IRAS beam-sizes overlaid, in 
Fig-El The IRAS measurements sample the complete detached 
shell, whereas the SWS spectrograph detects only a fraction of 
the shell. The different regions probed by the various instru- 
ments give rise to differences in the flux levels measured by 
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Fig. 7 . The detached-shel l model compared to the observations of R Scl in bl ack. The optical p hotometry is shown in trian- 
gles (jjohnson et alJ ll966l). Between 2 and 45 pm we show the SWS spectrum faonv et akl EoOZ): between 45 a nd 200 urn we 
show the LWS spectrum. The diamonds represent the IRAS measurements and the submillimetre points from Bagnul o" et alJ 
( 1998 ). The dashed gray l ine is the model SED, taking identical parameters for the dust shell as derived for the CO shell by 
Schoi er & Olofssonl bOOlh . The full gray line shows the model spectrum after correction for the limited size of the apertures of 
the SWS spectrograph. The dotted line represents the photospheric continuum. We also show the synthetic IRAS fluxes as derived 
from the model SED. The 25 and 60 pm IRAS fluxes are in excellent agreement. The predicted flux level at longer wavelengths 
is slightly too low. The model explains the observations very well. In the inset, we show the fit to the "30" pm feature on a linear 
scale. We also show the IRAS/LRS spectrum, which starts to deviate from the SWS spectrum at the wavelength, where the model 
predicts a difference, due to the different beam-size of IRAS and SWS. 



these instruments. The IRAS measurements sample the com- 
plete shell and should therefore yield much higher flux level 
than SWS does. This is indeed as observed (see Fig. [7). A 
smaller difference in the flux level, due to the size of the SWS 
apertures, is also observed in the SWS data below 27.5 pm and 
above 29 pm. 

The model is able to explain the strength of the MgS 
resonance within the SWS aperture accurately. It simultane- 
ously explains the much stronger excess in the IRAS beam. 
The presence of the MgS excess was already suspected by 
Bagnulo et al.1 dl998l) who found that the IRAS 25 pm pho- 
tometry was higher than expected on the basis of the flux level 
of their Cooled Grating Spectrograph (CGS3) spectrum which 
goes out to 23.5 pm. These authors speculate that the high 
IRAS 25 pm flux level might be due to a strong emission fea- 
ture, at wavelengths longer than 24 /mi, i.e. the "30" pm fea- 
ture. When we convolve the synthetic SED with the IRAS filter 



profile, the IRAS 25 pm flux is indeed nicely reproduced by 
our model. 



The excellent agreement, between the simple detached- 
shell model and the ISO and IRAS data, demonstrates that the 
"30" /im feature arises from the location of the detached shell. 
We conclude, that MgS grains are present in the detached shell 
of R Scl. We would like to emphasise that we do not attempted 
to further optimise the obtained model by varying the inner or 
outer edge of the dust shell. We have to make a rather crude as- 
sumption for the absorption level of the MgS grains below 1 /mi 
(Sect. l2~2l . Therefore, this simple model does not warrant such 
a detailed fitting procedure. However, the correspondence be- 
tween the model and observations, especially those probing dif- 
ferent spatial scales, demonstrates that the MgS emission arises 
from roughly the same location as the molecular gas emission. 
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Fig. 9. The single detached-shell model applied to U Cam. Symbols are the same as in Fig.0 We use the inner and outer radius 
and the gas mass of the CO shell and a 4% MgS fraction. The model that fits the SWS observations best does not reproduce 
the IRAS flux levels. Clearly, the single dust shell model is not applicable and the IRAS excess is predominantly due to more 
extended emission than the CO shell. 



4.2. Application to U Camelopardalis 

We also apply the model to U Cam. Compared to R Scl, this 
object has a much thinner CO shell (Table Q. Linda vist et al. 
(1999) measured an inner radius of 5.5 10 cm with a thickness 
of only 1.1 10 cm at a distance of 500 pc. This corresponds 
to an outer radius of 9" and the shell fits entirely within the 
largest SWS aperture. With those values for the inner and outer 
radius of the shell the model cannot explain the IRAS and the 
ISO measurements simultaneously (Fig. |5Jl. We find that the 
low IR flux levels, as measured by ISO/SWS, place a stringent 
constraint on the amount of dust located at the position of the 
CO shell. However, the IRAS fluxes clearly indicate that there 
is more dust present in the system. This dust is probably located 
even further away from the star. As can be seen in Fig. [5] all 
the IRAS fluxes as well as the 800 /mi flux are under-predicted 
by the shell model. When we increase the shell dust mass, in 
order to agree with the IRAS 25 and 60 //m flux levels, the 
longer wavelength fluxes are still very much under predicted 
by the relatively warm dust in the model. This indicates that the 
system contains cooler dust, which is probably located further 
away. The latter conclusion is in qualitative agreeme nt with the 
fact t hat the interferometric measurements of iLindqvist e t al. 
(1999) retrieved about half of the total CO emission, measured 



with single dish observations jNeri et alj |l998). This is due to 
a more extended component of CO emission. 

The situation in the case of U Cam is clearly more complex 
than for R Scl. The observations can not be explained by a sin- 
gle dust shell. The far-IR fluxes levels are dominated by dust 
located outside the detected CO shell. However, we stress that 
there is MgS present in the CO shell. The ISO/SWS instrument, 
which only probes the inner regions, detects a "30" //m feature. 
The temperature of the MgS, as derived from the centroid po- 
sition of the "30" yum feature in the ISO observation, yields a 
distance for the dust shell of 5.5 10 16 cm (7.3" at a distance of 
500 pc). This distance is identical to the size of the detached CO 
sh ell. Using the gas mass of t he shell of 1.4 1O~ 3 M as given 
bv lSchoier & OlofssorJfcOOll) and a dust mass of 2.4 10~ 6 M Q 
from the best fitting model for the ISO/SWS data, we find a 
gas-to-dust mass ratio of ^^SSO. 

We would like to point out that, in the two cases that we 
study here, the derived gas-to-dust mass ratio in the shell is rel- 
atively high (450 for R Scl and 580 for U Cam). These numbers 
are much higher than the values between 200 and 250, which 
are commonly used for modelling carbon-rich A GB and post- 
AGB stars (e.g. Ijurall 1 986t iMeixner et alJll997t) . However, the 
numbers we d erive are in relatively good agreement with the 
high values of Schoier & Olofssonl J2001I) . One should clearly 
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Source 


niboi 1 




R-skell ["] 


W Ori 


2.87 


31.3 


8 


SCep 


3.15 f 


31.2 


6 


TDra 


4.94 


30.2 


1.6 


VCyg 


3.24* 


30.5 


4 


RU Vir 


4.24 f 


30.4 


2 



Table 2. Estimated detached-shell radii (R S h e n) of warm car- 
bon stars exhibiting a cool MgS feature in their SWS spectra, 
based on the relation in Eq. between the "30" pm feature 
centroid position and th e size of the shell. 'Bolometric magni- 
tude from Bergeat et al. (2002). 2 Centroid positions taken from 
iHonv et alJJ2002h . Average of the listed values. 



be cautious in applying a "standard" value, when determining 
the total shell masses for such systems based on the IR excess 
alone. 

4.3. Other warm carbon stars with the "30" /urn feature 

There are a few more warm and optically bright carbon stars 
that exhibit a cool MgS feature in their SWS spectra. We do 
not model these stars here, mostly because of the lack of suit- 
able molecular data to compare with. However, we do give a 
crude estimate of their shell radii. We calculate the angular ra- 
dius of the star times the effecti ve temperature square d from the 
bolometric magnitude given bv lBergeat et alJ J2002h - using: 



Wboi,* = Mboi,* - 2.5 logi 



^★,arcsec T * 



Mboi,© - 5 log 10 
(A^boi,© - m bo i,*) 



7io 

Re 

10 



2\ 

j 

To 



M b , 



•Ol,0 



— ) 

djj 



AfboLc 



10 AU 



:7V 



(3) 



where Wboi and Mboi are the apparent and absolute bolometric 
magnitude and /?*, arC sec is the angular radius of the star. The 
measured values of the apparent bolometric magnitude and the 
centroid position of the MgS feature together with Eq. ifj] & 
[3} yields the angular radius of the shell. The sources and the 
estimated radii are listed in Table [5] These sources all have 
much weaker "30" //m features than the feature observed in the 
SWS spectrum of R Scl, with a correspondingly lower signal- 
to-noise ratio. At the same time, the estimated stellar contin- 
uum underlying the "30" pm feature becomes more influential 
on the derived feature profile. Based on this, we estimate an 
uncertainty in /lc,30 of >1 pm and in the derived shell radii of 
at least a factor 2. 

5. Conclusions & future prospects 

With the models that we have presented, we now have a tool 
to systematically investigate the influence of MgS on the IR 
spectra of carbon stars with detached shells. We find that, in 
the case of R Scl, the dust (containing MgS) is co-spatial with 
the molecular gas. This shows that the dust present in the shell 
is indicative of C-rich chemistry and therefore the photosphere 



was C-rich at the time the shell was ejected. This rules against 
a scenario in which the shell is produced due to the transition 
from O-rich to C-rich. We find that the shell around U Cam 
contains MgS as well, again demonstrating the C-based chem- 
istry of the shell. 

The prominence of the "30" pm feature and the broad 
wavelength range over which the MgS grains emit make this 
material an excellent tool for future studies of the circumstel- 
lar shells of carbon-rich evolved stars. In particular, we point 
out that, because of the temperature dependence of the MgS 
emission profile, some of the usual ambiguity in interpreting 
SEDs can be lifted, as we can constrain the location of the shell 
by means of the "30" pm feature. Th e low-resolution lq ng- 
wavelength m ode of the IRS instrume nt ( Roellig et alJ 19981) on 
board SIRTF (Gallagher£tiD[2003) will be sensitive enough 
to easily detect the MgS feature in carbon stars located in the 
SMC or LMC. This will eliminate the distance uncertainty, 
which usually troubles the derivation of the shell parameters 
of galactic carbon stars and therefore allows a more statistical 
approach for studying the occurrence rate and shell masses of 
such detached shells. 

Another tantalising prospective is offered by the 
FORECAST instrument dKeller et alJ l2000l) for the SOFIA 
observatory (Beckl]3 119971) . With this camera it will be 
feasible to directly image the shells at the wavelengths where 
the MgS is emitting. The shells discussed here are all easily 
resolved with the good angular resolution (< 5") of the large 
SOFIA mirror. The several bandpasses in the "30" pm feature 
will allow us to study the MgS abundance and temperature as 
a function of distance to the carbon star to an unprecedented 
level of detail. 
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